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Magnetotransport with two occupied subbands in a Si(100) inversion layer 
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We have studied an electron transport in inversion layers of high-mobility Si(100) samples. At high 
electron concentrations and temperatures below 4.2 K, two series of Shubnikov-de Haas oscillations 
have been observed. The temperature damping of the second series oscillations indicates that 
the second occupied subband belongs to the first energy level of the fourfold-degenerate ladder 
0'. Samples with two occupied subbans exhibit a strong anomalous negative magnetoresitance, 
reaching « 25% of a zero field value at B = 12 T. The resistance decrease is more pronounced for 
lower temperatures and higher electron concentrations. We explain this behaviour by an increase of 
the second subband mobility due to the freezing-out of the scattering of 0' electrons. Based on the 
measured periods of SdH oscillations, we conclude that the electrons are distributed inhomogeneously 
beneath the sample gate. 
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The rt-type inversion layer on a (100) surface of p-type 
Si has been the subject of investigation for many years 
(for review see §). The silicon is an indirect gap semi- 
conductor having six minima (valleys) in the conduction 
band along the (100) direction. In the case of (100)- 
oriented wafers, two of the valleys in the Brillouin zone 
are oriented perpendicular to the Si/Si02 interface (with 
m z = 0.916), while the other four are parallel and have 
m z = 0.19. A quantum well is formed at the interface by 
applying a positive voltage V g to the gate, and the elec- 
tron motion in z direction is quantized. The difference 
in m z causes a splitting of the energy spectrum into two 
independent ladders of levels: a twofold-degenerate lad- 
der denoted 0, and a fourfold-degenerate ladder 0'. Due 
to the higher effective mass in the z direction, the lowest 
energy level in the potential well belongs to the 0-ladder. 

Disagreement exists in the literature on the electron 
concentration at which the Fermi energy crosses the next 
energy level and on the question whether it forms the 
bottom of the second subband of the 'unprimed' lad- 
der, or the first subband of the 'primed' ladder pl[||. 
More recently this problem was discussed in M, based 
on the combination of measurements of the gate voltage 
dependence of the magnetoresistance and the Hall volt- 
age. It was shown that the mobility of the second sub- 
band is higher than in the first subband (but relatively 
low), and that it linearly decreases with the tempera- 
ture. This temperature dependence was tentatively ex- 
plained by the temperature-dependent screening of the 
elastic scattering in the second subband. It was con- 
cluded from the experimental data that the second sub- 
band is the lowest subband of the fourfold-degenerate 
ladder, but the corresponding Shubnikov-de Haas (SdH) 
oscillations have never been observed explicitly. Here we 
report, to our knowledge for the first time, the magne- 
toresistance measurements on high-mobility Si(100) sam- 
ples, showing clearly both series of Shubnikov-de Haas 



oscillations and supporting the idea that the second oc- 
cupied subband belongs to the first energy level of the 
fourfold-degenerate ladder. 

We employed in our experiments five Hall bar samples 
prepared from three different wafers. The wafers were 
made from p-type silicon with resistivity 20f2cm at the 
room temperature, their crystallografic orientation was 
checked by X-ray diffraction and found (100) with an ac- 
curacy within ±0.5%. The edges of samples, defined by 
the gate geometry, were oriented along (110) direction. 
Results of measurements were very similar for all sam- 
ples used. Therefore, we present data obtained for two 
of them, denoted as 'sample 1' and 'sample 2'. The max- 
imum mobility of electrons in 2D inversion layers was 
about 2.7 m 2 /Vs for the sample 1, and 2.0 m 2 /Vs for the 
sample 2 at the temperature T = 4.2 K. The dimensions 
of samples were 5 x 0.8 mm 2 , with the distance between 
potential leads 2.5 mm, (sample 1) and 2.5 x 0.25 mm, 
with the distance between leads 0.625mm (sample 2). 
The 0.2 /im thick SiC>2 gate insulator allowed to feed 
the electrons into an inversion layer with the maximum 
concentration N s as 1.3 x 10 17 m -2 . There is a linear de- 
pendence between AT, and V g determined by the sample 
capacitance; the rate of filling, N s /V g , was approximately 
1.1 x 10 15 m _2 V _1 for both samples. Shubnikov de Haas 
oscillations were investigated either as a function of V g , 
for a constant magnetic field B, or as a function of B, 
assuming fixed V g . We used a DC measuring method, 
the current through the sample was varied in the range 
0.1-100/kA, and the voltage between the potential leads 
was measured by a digital voltmeter. The samples, im- 
mersed in the liquid 4 He, were cooled down to 1.9 K in 
magnetic fields up to 12 T. 

Figure |l| presents the selected measurements of mag- 
netoresistance as a function of the gate voltage V g . 

Firstly, the upper panel shows the 'standard' data ob- 
tained using the measuring current / = 10/iA; these re- 
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suits agree quite well with results of previous measure- 
ments g . The observed oscillations of magnetoresistance 
correspond to filling of Landau levels of the first subband 
and they are periodically dependent on V g for the electron 
concentration lower than w 8 x 10 16 m~ 2 . Due to the spin 
degeneracy, and the twofold valley degeneracy, there are 
A\e\B/h electrons per unit area in each Landau level. The 
amplitude of oscillations drops down rapidly for higher 
electron concentrations, and their period decreases. To 
help resolve the oscillations, the inset presents the deriva- 
tive of the magnetoresistance curve with respect to V g for 
higher N s region. 

10" 16 xN s (m~ 2 ) 

2 4 6 8 10 

2.5 i ■ ' ■ ' ■ ' ■ ■ ■ ' 1 



^2.0 




10 20 30 40 50 60 70 80 90 100 

V e (V) 



FIG. 1. Shubnikov - de Haas oscillations and filling factors 
for sample 2 at B = 4.18 T, I = 10 fj,A, and T = 1.9 K. The 
upper right curve is a derivative of a main one with respect 
to V g , with an arbitrary j/-scale (thin lines). The thick line 
shows results of a measurement using a reduced current 0.1 
fiA. 

The lower panel shows the filling factor of electrons 
in the first subband, determined from the minima of 
the SdH oscillations shown above. For V g < 80 V, the 
concentration N s> \ follows the straight line, implying 
N St \ = N s . For higher V gi a single period oscillations re- 
main present, but the electron concentration in the first 
subband is lower, since a part of electrons is fed into the 



higher subband. 

Secondly, a part of an 'anomalous' magnetoresistance 
curve is presented in the upper panel of figure @. The 
measuring procedure was the same as described above, 
but the current through the sample, /, was reduced one 
hundred times, to 0.1 /iA. In that case the second series 
oscillations became visible and the decrease of the mag- 
netoresistance is observed. It is important to stress that 
the resistance measured at B — was identical for both 
measuring currents, i.e. the same curves were obtained 
for / = 10 /jA and I = 0.1 /iA. Therefore, it is not pos- 
sible to explain the drop of magnetoresistance simply as 
due to the reduced heating of 2D electron gas. If it would 
be the case, the reduced electron temperature should im- 
ply e.g an increase of amplitude of SdH oscillations of the 
first subband electrons. Instead of it, the decrease of am- 
plitudes was observed, as if a part of a current, originally 
carried by the first subband electrons, was taken over by 
the second subband electrons with higher mobility. This 
indicates that only the electrons from the second subband 
are influenced by the the current magnitude. 
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FIG. 2. Magnetoresistance oscillations for several electron 
concentrations and two temperatures. Sample 1, / = 5 /J.A. 

These unusual results motivated us to more labori- 
ous investigation of the field dependence of magnetore- 
sistance with the fixed gate voltage V g and the variable 
magnetic field. In that case the coexistence of two series 
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of oscillations can be studied, together with the devel- 
opment of the negative magnetoresistance. Examples of 
such curves, obtained for four gate voltages at two differ- 
ent temperatures, and employing a relatively high current 
1 = 5 /iA, are presented in figure ^. In the lowest mag- 
netic fields, the negative magnetoresistance, attributed 
to the localization, is observed for all four V g , followed 
by a slight increase, which lasts till B « 3 T. In this low 
field region, the magnetoresistance practically does not 
depend on the sample temperature. For higher fields, 
the SdH oscillations appear, corresponding to both sub- 
bands, and their emergence is accompanied by onset of 
a strong negative magnetoresistance. The field depen- 
dence of smooth parts of curves is very sensitive to the 
temperature, for lower temperatures a more pronounced 
resistance decrease occurs, similarly as in the case of re- 
duced current density, as presented in figure |l|. 




are shown in figure ||. The lines are offset for clar- 
ity. Peaks corresponding to two different frequencies are 
clearly seen. The first subband oscillations are almost si- 
nusoidal and, therefore, described by a single peak which 
has somewhat lower amplitude for the curve measured 
at lower temperature. The second subband oscillations 
contains a large number of higher harmonics and the cor- 
responding peaks are more visible at lower temperature 
curve. This can be understood having in mind a semi- 
elliptic form of oscillations. The smooth background is 
due to the smooth part of the negative magnetoresistance 
which was not removed before transforming the data. 
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FIG. 4. Electron concentrations calculated from periods of 
magnetoresistance oscillations. Also a difference between the 
full concentration N s , and the first subband electrons concen- 
tration jV a ,i is shown. Sample 1, I = 0.1 jj,A. 



FIG. 3. Absolute values of Fourier transforms of high-field 
parts of magnetoresistance curves measured at V g = 99 V and 
shown in figure |^. 

The magnetoresistance curves, measured as a function 
of B at various temperatures, allow us to derive im- 
portant characteristics of the second subband. Firstly, 
we have determine the cyclotron effective mass from the 
temperature damping of the oscillations amplitudes, and 
found m c = 0.43 ± 0.03, in good agreement with a the- 
oretical value, m c = 0.417. This implies that the first 
subband of the 0' ladder is occupied, in agreement with 
conclusions of references and || . 

Both series of oscillations are periodical functions of 
the inverse magnetic field and their frequencies are lin- 
ear functions of concentrations of carriers in the individ- 
ual subbands. Examples of Fourier transforms of two 
curves, measured at V g = 99 V and presented in figure 



Both N a ,i and N Si 2 determined from the periods of 
SdH oscillations are shown in figure There are only a 
few points of N Sy2 (V g ) dependence and, therefore, the lin- 
ear extrapolation to lover gate voltage has only a limited 
meaning. Nevertheless, it indicates a lower threshold of 
occupation of the second subband than reported previ- 
ously. The results also show that the ratio of two periods 
of oscillations is approximately 4.5 for all four investi- 
gated gate voltages, which means N St i/N s> 2 ~ 2.25, if the 
fourfold degeneracy of the second subband is assumed. 
Traditionally, ||-[|] , the concentration of electrons in the 
second subband, N s> 2, is deduced from the decrease of 
the period of the first subband oscillations at higher gate 
voltages, such as shown in figure |l|. If the homogeneous 
distribution of electrons beneath the sample gate is as- 
sumed, iV Sj 2 is given by the difference N s — N s ^. Our 
direct determination of N Si i and N s ^ from the oscilla- 
tion periods does not confirm results of this approach in 
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the case of N Si2 - 

More series of SdH oscillations were already reported 
for inversion layers on surfaces slightly tilted from (100) 
H and for layers on (100) surfaces subject to cither uni- 
axial |6) or biaxial stress 0. In the first case, only the 
first subband electrons are involved and the reason for 
additional periods is lifting of the twofold valley degener- 
acy caused by the tilting and the inter-valley interaction. 
In the latter case, the crystal symmetry is lowered due 
to the stress and, consequently, the energy of 0' levels 
shifted below the Fermi energy level. 

We do not think that our second series oscillations can 
be satisfactory explained by the above mentioned mech- 
anisms. The surface orientation is (100) with high ac- 
curacy and the only sources of stress are the different 
thermal expansion coefficients of Si, SiC>2 and the gate 
material. On the other hand, we would like to turn at- 
tention to another possibility, not mentioned previously, 
why the bottom of the first subband of the 0' ladder can 
be shifted below the Fermi energy only in a part of the 
sample area: the inhomogeneity of distribution of elec- 
trons beneath the sample gate. 

In our samples, the area occupied by a two-dimensional 
electron gas is defined by the gate geometry rather then 
by fixed ionized donors, as it is the case in modulation 
doped heterostructures. The concentration of electrons 
N s is related to V g by N s = eV g /\e\d where e is a di- 
electric constant of the oxide and d is its thickness. The 
charged gate and the 2D electron system act as a paral- 
lel plate capacitor; the plates containing the free charge 
keep the constant potentials through their areas, rather 
then the constant concentration of carriers. Therefore, 
the free charge density can fluctuate inside the sample, 
e.g. due to the surface roughness and/or an inhomoge- 
neous distribution of charged impurities in the depletion 
layer. Moreover, the free charge should accumulate near 
the plate edges, to compensate the Coulomb repulsion of 
carriers from the sample interior. The electrostatic cor- 
rections to the finite dimension of a capacitor [0 predict 
that the extra charge is localized in narrow strips near 
the edges, with the width proportional to d, and that 
the excess concentration, AN s (x), diverges when x, the 
distance measured from the edge, approaches zero. In 
real life the concentration will stay finite, nevertheless 
we can assume that it first rises near an edge, before it 
vanishes. Thus, we consider that the narrow channels 
near the edges concentrate a part of electrons from the 
second subband and contribute substantially to the sam- 
ple conductivity. 

The high sensitivity of the observed results to the mea- 
suring current density and to the temperature suggests 
that the electron-electron interaction can be involved. 
Firstly, the temperature-dependent screening of the sec- 
ond subband electron scattering can play a role, as de- 
scribed in 0. Secondly, the properties of the second 
subband electrons should be strongly influenced by quan- 
tization into the Landau levels for hu> > ksT ' . Note that 
the onset of the negative-dependent magnetoresistance 



occurs just when the above condition is satisfied. 

The Landau levels of a single subband never cross and 
their shape more or less copies a shape of a smooth confin- 
ing potential and exhibits a characteristic step- like struc- 
ture near abrupt potential steps. The Landau levels of 
electrons from two different subbands, and 0', can cross, 
as they will see the shape of a sample potential in a dif- 
ferent way. Firstly, due to a large difference between the 
number of Landau levels near the Fermi energy and, sec- 
ondly, due to the difference in the effective masses. The 
points of crossing can represent the 'hot' points of an in- 
tensive intersubband scattering which can substantially 
contribute to the observed effects and can be suppressed 
at low temperatures and high magnetic fields. 

The increase of the mobility of the second subband 
electrons causes the redistribution of the current den- 
sity and more current is conducted by them. As a 
consequence, the amplitudes of second series oscillations 
growth to detriment of amplitudes of the first series os- 
cillations. The extension of present measurements down 
to the milikelvin temperatures region, which is now in 
progress, can probably shed more light to the nature of 
this new interesting effects. 
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